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INTRODUCTION 

Recent e>:perirnenta1 and theoretical studies related to ultrasonics 
have demonstrated that the scattering of elastic waves by defects in a 
material provides information about both. This information may be used 
for determining the properties of any imperfection such as the size and 
orientation of a crack or it may be used for determining the mechanical 
and strength properties of the material. 

Vary and his associates have emphasized the aspects of applying 
ultrasonics to material evaluation [1,2,3]. Their material information 
contents are experimentally measured in terms of ultrasonic attenuation 
and velocity factors. Using these concepts and a simple model, Vary [4] 
presented some very useful empirical relations that correlate these factors 
to the fracture toughness and the yield stress and gave data for two 
maraging steels and a titanium alloy. 

The above mentioned findings are indeed very interesting and yet 
somewhat intriguing. After going through an extensive literature review, 
the author [5] suspected that the links between the strength properties 
(fracture toughness and yield stress) and the ultrasonic factors (attenu- 
ation and velocity factors) are the material microstructural parameters. 

It was pointed out that the size of the second phase particles and the 
distance between them play a very important role in material resistance 
to fracture. 

In a first step toward defining the relation between , fracture 
toughness and Vjg^/m , the ultrasonic attenuation factor, the interaction 
of a pair of second phase particles (referred to as inhomogeneities) is 
studied. Since energy is trapped or dissipated in the vicinity of the 
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neighboring inhomogeneittes , the ultrasonic factors measured at far field 
reflect this energy loss. Assuming that the attenuation is a function of 
incident wave frequency and that the energy trapped in the vicinity of the 
inhomogeneities provides the fracture energy a relation between the fracture 
toughness and the ultrasonic attenuation factors can be obtained. To ac- 
complish this purpose both the far field and the near field solutions of 
the interaction problem are needed. 

Currently, scattering theory of a single flaw is available [5], The 
approach of Gubernatis, Domany and Krumhansl [6] is of particular interest 
in that they gave the scattered field solution far from the flaw in terms 
of the displacements and strains in the scatterer, i.e. the inhomogeneity. 
They also expressed the physical quantities such as total and differential 
cross sections in terms of the scattered field quantities. Since the 
strains and di spl acements inside the scatterer are not available for the 
dynamic case, they used the results from Eshelby [7] for the static case 
and studied the case appropriate to long wave limit [8]. 

The purpose of this report is to study the dynamic response of in- 
homogeneities, one or two, and to determine the strains and displacements 
inside them under incident plane waves as depicted in Figs. 1 and 2. The 
method of equivalent inclusion is used [9]. The underlining approach, 
the formulation and governing equation for the ei genstrai ns , and the 
determination of the energy due to the presence of the inhomogeneities 
are presented in this report. The derivation of the correlation between 
K,, and v. 8 will follow later. 

It. L ,5 
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RE VIEW OF SCATTERI NG THEO RY: SINGLE FLAW 

In a recent article, Gubernatts, Domany and Krumhansi (GDK) [6] 
summarized their recent work on elastic wave scattering with application 
to nondestructible evaluation. They studied the elastic scattering of 
a single flaw Che it a void, a crack or an inhomogeneity) under an experi- 
mental situation as depicted in Fig. 1. The scattered amplitudes and 
cross sections (measurable quantities) were derived specifically in terms 
of the scattered fields of stress and displacement at large distance from 
the scatterer. 

For an incident plane wave of angular frequency co and wave vector IT , 
Fig. 1, the incident displacement vector is 

u^.(r.t) = u^ exp i(k • r - cot) (1) 

where 

|kl = a = co/V|_ 

B = (o/v .j. 

V|_ = longitudinal sound velocity 

Vj = transverse sound velocity 


Employing integral theorems and the Green's function approach, they 
showed that the final form of the basic scattering equation is 


UjCD = u,-(r) + U|(r) 


( 2 ) 


in which 


u^.(r) = $p co^ / dv‘ g^-^(r - r’)Ujj,(r') 


where 


( 3 ) 



4 


’^flaw Matrix 

^ " ^‘^jklm^flaw “ ^^jklm^matrix 

^2 ■ region occupied by scatter 


For an isotropic elastic medium, in the far field at a distance r 
from the defect, they found that the scattered field depended on a certain 
vector, the f-vector, as follows: 


^lar 

“i ■ “ > F + 


* t«,j - =,=j) fj( 6 ) ' 

where the f-vector is defined as 


igr 


(4) 


f-(k) = — [5p 0 )^ / dv u- exp(-ik*r) 

1 AttHM rj 1 ‘ 


Attpo) 

+ i k e. 


j ® S'jkl ^ ^kl exp(-ik-r)] 

R, 


(5) 


and is dependent upon the differences in material density and properties, 
between the matrix and the defect, the incident wave field and the total 
displacement and strain fields "inside" the scatterer, R 2 . 

With the asymptotic value of u^ given in Eq. (4), the asymptotic 
value of a|j can be obtained as 

piar 

o. . ~ iXa - e^ 6. . f^^(^ 

eiar _ 

+ iy [2a - [e.e^ej^f|^(a)] 

pi^r _ _ 

+ [e.f^(B) + e^f.(3)] 

piPr _ 

- 23 - e.e.e f (3) ) r -*■ “ (6) 

r 1 j k k' 
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and various cross section can be expressed in terms of the scattered 
asymptotic values. For example, the differential cross section, for any 
given frequency it is a measure of the fraction of incident power 
scattered into a particular direction, is found to be 


dP((iJ) _ lim <r^e-j a?- uj> 
dfl <I°> 


(7) 


where a dot denotes time differentiation and an angle bracket denotes 
time averaging, i.e. 

<f(t)> = Y • f(t)dt 

0 

The differential dfl is the differential element of a solid angle. The 
total cross section is simply 

P(ca) = / dn (8) 

4tt 

From Eqs. (4-8) it can be seen that the far field solution, the 
scattered amplitudes and cross sections, depends upon the determination 
of the displacement and strain fields inside the scatterer and the evalu- 
ation of the volume integrals in Eq. (5), see statements following Eq. 5, 
Ref. [9]. The displacement and strain fields are currently not available 
for the dynamic case. As a remedy, Gubernatis [8] used the results 
obtained for the static case from Eshelby [7], obtained by employing a 
method called equivalent inclusion method. 
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INHOMOGENEITIES IN A TIME-HARMONIC HAVE FIELD 

From the above formulation, it is clearly seen that displacement and 
strain fields inside scatterers (inhomogeneities) play a very important 
role in nondestructive evaluation. In studying the interaction between 
two inhomogeneities, the dynamic version of the equivalent inclusion 
method is used. The method was first employed by Mura [10] and his 
associates 111] in studying composites. 

Let the ultrasonic experimental situation be depicted as in Fig. 2 
such that the incident power is along the positive z-axis. Let the total 
strain be the combination of elastic strain and non-elastic (or eigen-) 
strain : 


e 

rs 



(9) 


P ★ 

where and are the total strain, elastic strain and the 

eigenstrain, respectively. The eigenstrain is the non-elastic strain 
which is caused by a change of the form of an inclusion, which if the 
surrounding matrix material were absent, would have gone through some 
homogeneous deformation. Due to the difference in material properties 
between the matrix and the inclusions, this change of form causes dis- 
turbance in stress and associated strain in the material. 

The equations of motion for a continuum are 


ajk,k = P Uj (10) 

where a dot indicates a differentiation with respect to time while a 
subscript comma indicates a spatial differentiation. In a linear elastic 
matrix with small strain deformation, the strain-displacement relations 


are 
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rs 


2 ‘^s.r^ 


(.11) 


and the generalized Hooke's law is 


®jk ^jkrs ^rs 


( 12 ) 


where are the elastic constants. Employing Eqs. (9-12) the dis- 

placement equations of motion can be written as: 


^jkrs “r.sk ” ^jkrs rs,k 


= c 


(13) 


If the matrix is isotropic the elastic constants can be expressed 


as 


^jkrs ^ '^jk'^rs ^ '^jr'^ks ^js'^kr 


(14) 


where X, y are Lame's constants and 6^-j is the Kronecker's delta. 

If the displacement and strain fields are time-harmonic, they can 
be written as 

(15) 


Uj(r,t) = Uj(7) exp(-iojt) 


e^ 5 (r,t) = er*(^ exp(-iwt) 


(16) 


where oj is the frequency of the incident wave and i • i = -1 . A substi- 
tution of Eqs. (15,16) in Eq. (13) leads to 


'^jkrs ^,sk + = ‘^jkrs %s.k '' 


(17) 


For a body of volume v and surface s containing an arbitrary distribution 
of eigenstrain, the traction free boundary condition can be expressed as 


^jkrs ‘^CjS "k 


‘^jkrs ^rs "k 0" " 


(18) 
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\ 


Consider now the associated Green's function with homogeneous 
boundary conditions: 


<=jkrs 9rm.sk - P 9jm = ’ ^jm " 

^jkrs \m,s \ = 0 , on s 


09 ) 

( 20 ) 


Here, the Green's function gjjjj(r-r ' ,t) represents the displacement in the 
j-di recti on at point r by a unit body force in the m-di recti on applied at 
the point r' . The points defined by r and F' are referred to the obser- 
vation and source points, respectively. The Dirac delta function, 
6(r-r',t), represents the body force. For the time-harmonic case 


6(r'-r',t) = 6(r-r') exp{-io)t) (21) 

gjn,(r-F'.t) = gj^(F-F’) exp(-icot) (22) 

The substitution of the above equation in Eqs. (19,20) yields 

^jkrs 9rm,sk P"-'" 9jm = " jm (23) 

"jkrs S^m.s "k = ° ^ (24) 


To derive an integral representation for the solution to the eigen- 

strain problem, a dynamic version of the Betti-Rayleigh reciprocal theorem 

is considered. Multiplying Eq. (17) by g. and Eq. (23) by u. and sub- 

J*” J 

stituting one from the other the following is obtained: 

'^m(^ ^ ^jkrs (9jm “r,sk " “j 9pm, sk ~ 9jm ®rs,k^‘^'' (29) 

after an integration over the volume v. Noting the symmetric properties 
of and applying Gauss' divergence theorem and boundary conditions, 

an integral expression of the displacement field in terms of the eigen- 
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strains is obtained as follows; 

*'jkrs ^jm'^rs.k 

With the application of Gauss*^ theorem once more, the above equation can 
be reduced to a more compact form: 


u^tP ) 

‘^jkrs ^jm,k 

(r-r') e^*(F)dv 

(27a) 

or 

u„,CF) 

m 

■J ^jkrs 3jm,k 

Cr-r’ ) CpsCr' )dv' 

(27b) 

where 


9 

^jm.k ~ 9X|^ 9j[j, 

_ _ 9 

9xj^ ^jm 

(28) 
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ISOTROPIC LI NEAR ELASTIC MEDIUM 

For a linear elastic, isotropic and infinitely extended material, 
the Green's function is well-known: 



(r-r') 


__L_ r .2 exp iBR 

47Tpo)^ R 


where 


8 8 , exp iaR exp ipR , 

8x . 8x ^ R R ^ ^ 

J ^ 


R = |r-r'| 


C29) 



By substituting Eq. (29) in Eq. (27b) the induced displacement can 


shown to be 


u (r) = - j [ a^X Y 

m 4TTpo) rr,m 


^ 'J'jk.jkm 


where 


* ^Jk.jlon ‘ 


- I" 

Q 

♦u -m 


e,*( 7' ) HPOaEl dp 
E.J(P) «P(Mi dP 


(30) 

(31) 
(3?) 


From Eq. (30) the strain field can be obtained by direct differentiation 
as follows: 


J 



n 




- T--T 

4Trpo) rr,mn 

^ ^ ^k,km^ 

'*^jk,kjmn '*' *jk,kjmn ^ 


(33) 


The displacement and strain field given in Eqs. (30,33) are those 
due to the presence of eigenstrains in the regions of inhomogeneities 
n's, respectively. They obviously depend on the form of the eigenstrains. 
Since the eigenstrains are not a known priori, it is convenient to expand 
the eigenstrain in a form of polynomial 1 12] 


E.. - B.. + B..| X, + B..|, X. X, + ... 

ij ij ijk k ijkl k 1 


(34) 


in the region where the eigenstrains are present. The quantities B.jj, 

B.jjl^, ... are constants symmetric with respect to the free indices i and 
j and having values independent of the order in which the summation indices 
appear, i.e. , etc. Using Eq. (34) for 

e.jj and substituting it in Eqs. (31,32), the functions T.jj and O.jj are 
found to be in terms of the constants B.jj, ... and some volume 

integrals as follows: 


4'.. = B-. ^(r) + B..^ ^^(r) + B.j^., 4^^^(r) + ... (35) 

= B.. <(>(?) + B.j^ 4>k(r) + d>^^(7) + ... (36) 

where 

^(r) = /f/ ^xP( . i^R ) dv' 

'Pk('r) = //J dv' 
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M...S 

♦(7) . JJJ dv' 

n 


C37) 

(cont'd) 


♦k(F) . Ji; X' dv 


kl 


(r) = HI x' X' ... 

, . .S O K 1 


, exp(igR) 


The substitution of Eqs. (35,36) in Eq. (33) leads to 


®kj °mnkjl^^ ®kjl ' 


(38) 


where B. . are constants, and 

Kj . . . S 


- 6kj) 

4iTpco Dj^n(^j-l(r) - 


(39) 


ug^(<t>T I 6. +(I)t 1 6.) 

^ ^^1 ,kn jm ^1 ,mk jn' 


^ 'i^l ,nin ®kj 


It should be noted that the are symmetric with respect to k,j and 

m,n. Generally, D-|^i_ _(r)?‘ (r) unless ij^j and k^l , e.g. 

° n 22^°2211 ° 1223^°2312 ’ 

The development so far has reduced the determination of the displace- 
ment and strain field to the determination of the constants __ 
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and the volume integrals given in Eq. (37). For the static case, 
a , g 0, Eshelby showed that Dninkj^^^ ® constant and Dninkjl . .. 
zero. Once the integrals in Eq. (37) are evaluated, the solution now will 
depend upon the determination of the constants 
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EVALUATION OF VOLUME INTEGRALS 4>‘S AN D ip'S 

Let the volume integrals given in Eq. (37J be denoted by 


I -m p(P) • 


(41) 


where C = a for the ijj-integrals and C = 3 for the <})-integral s , p(r') is 
of the form of (x')^ Cy')^ (z')^ , is an interior region where p(r') is 
distributed, and R = |r - r'|. Employing suitable Taylor series expansion 
and the multinomial theorem, the I-integral can be written in a reasonably 
convergent series as 

“ n n-1 /- Tvn „n 

_ Y Y Y (-1) 9 /-exp icr^ 

^ n=o 1=0 k=o llk:(n-l-k)l ‘ 9 xi 9 yk 9 zn-l-k ^ ^ ^ ‘ 


/// (x')^y')'^(z')"'''''' p(x’,y’,z')dx‘dy'dz' , r outside (42) 

Q 


and 


n n-1 (_i )'^ 


!(<")< - I I I l!k!(n-l-k)! * ^ 

n=o 1=0 k=o 


1 .,k^n-l-k 


9 " 






r inside Q 


(43) 


When is an ellipsoidal region the integrals in I^ and I^ can be 
readily obtained by using results from Dyson [12], e.g. 


Ill — dv' = 4ua a a J ^ ^ 

n ^ ^ m=l (2m-l ) ! (2m-l ) (2m+l ) 


f , ,m-l 2m-l m-1 m-1 

(-U c ai. aj£ 
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III i, ^'-Ssrk- • 

j. 12m-2)(2m-4) m-3 m-3, . (2m- 2)(3) 1 

f r2i-“3Tr25^ ’k ''p r25i-3)(2m-l) =k =k ^ 


III “ „i 

' m=o 


. y (-1)" 1 ^ 

2 m! fn,o 


where 


iraiajaa 

‘■'".0 ' PSTmTlT 


2mi 2m2 ^tUg ! 2mz ! ! 

® ^ ^ mj im2 img ! 


° (a, H*)"'’(a,'+*)”“(a, “+♦)"" 


my + m 2 + [TI 3 = m , 

Q = (aj^+i(;)(a2^+^ij)(a3HTlj) , 


in which (x/aj^ + (y/a^)^ + (z/a^)^ = 1 


1 
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DETERMINA TIO N OF EIGENS TRAINS 

Let the inhomogeneities be situated as shown in Fig. 2, where 

and fljj are the regions occupied by the inhomogeneities. If the 

applied elastic fields are denoted by strains c^j(r), stresses ^^^(r) 

and displacements u^(r) and if the self-equilibrated fields due to the 

presence of the inhomogeneities are denoted by e- -(r), o-.j(r) and u.(F) 

I j ' j j 

the total fields due to the applied field and the inhomogeneities are 

then the sum of the two. It is noted here that the time-dependence is 

•k I 

suppressed. The MEI can be used to determine the eigenstrains (r) 


in xQj and 

★ 

e. . 
ij 

II,-x . 

(r) in 

fill 

as foil 

ows : 




^^ijkl ^kl 

(F) 

''ijkl 

* 

^kl 

^(F) = - 

ACi^.ki e 

",l(F) i, 

1 flj 

(44) 

^^ijkl ^kl 

(F) 

0 

+ C • • 1 1 
1 jkl 

★ 

^kl 

^^(F) = 

-AC 

^S-jkl 

£ki(F) 

in fljj 

(45) 

where 











I 


I 

0 






AC. , 
ijkl 

“ 

^ijkl ■ 

^•jkl 



(46) 



II 


I 

0 



(47) 



^^ijkl 


^•jkl 

^i jkl 



0 


I 

II 






and C. , 
1 jkl 

, c 

ijkl ’ 

ijkl 

are the 

elastic 

moduli 

tensor of the 

matrix. 


the first inhomogeneity , and the second inhomogeneity Iljj, respectively. 

Consider now the expansion of the eigenstrains in polynomial form 
such that 


-k X, ♦ 

e*/( 7) • B,] t B,],^ + B.'^, X, t ... in a, 


( 48 ) 


( 49 ) 
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_* II,-, II II _ II _ _ 

E. . Cr) = B . . + B . X, + B . , X X, + . , . in Q.. 

1.1 IJ Ijk k ijkl k 1 II 


(50) 


where the barred and unbarred quantities are measured from the two 
different coordinate systems as shown in Fig, 2. The position vectors 
are related by 


^ = (x*^ - x°) a 

X. Uj- X.) a^. 


(52) 


where a., is the coordinate transformation matrix. It is easy to show 

I II I II , , 

that the unknown constants B— , B^j , , B^-jj^ , ... in Eqs. (49,50) 

satisfy the following simultaneous equations: 


^^stmn ^ ^^mni j ^®^^i j '^mni jk ^®^^i jk 

'^mnijkl f°^^ijkl ^mc^nh ^^chij^^^^ij 


II 


II 


II 


II 


" ^chijkf°]^ijk " °chijklt°^^-jkl " 


^stmn ^mn " "^^stmn^mn 

^^stmn ^^" 5 ^ ^mnij^^^^ij - 5 ^ ^mni jk ^®^®i jk 
^ “^mni jkl f^^^ijkl 


®mc ^nh ®pf °chijf°^^ij 

II 


II 


^ 3x^ '^chijk'^^^^ijk ^chi jkl ^®^®i jkl 

0 I _ I I 

^stmn “mnp “ ~^^stmn ^mnp 


]} 


(53) 


( 54 ) 
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21 ^^stmn 3x °mnl j j ^ 3x^9x '^mni jk 

r H r H 


^ 9Xp3Xq "^mnij-kl ^°^^ijkl + ••• ^ 


+ a a u a ^ a r— ? D u - ■ [ 0 1 B - - 

°mc nh “pf “qg ^3xf3x chij^ J ij 

I □ 


^ 9x^3Xg ^chijk^^^^ijk 

92 II II 

31 Tf 3 'Xg ^chi jkl jkl ''■ ••• 


0 I I 

+ r R “ AT F 

stmn mnpq " “ ^stmn ‘^mnpq ’ 


'^^stmn ^^cm%n '*' ^chi jk jk 

’^chi jkl jkl •••^ ■'■ t^mnijtOlBij 


^stmn ®mn “ 


stmn ‘-mn 


^^stmn ^^cm^hn^fp °chij^°^^ij 

'*' ‘^chi jk ^ ^chi jkl ^°^®ijkl 

II II „ II II 

+ [_^ D . .[0]B. . + _i_ D -.^[OlB.., 
^3)Tp mnij ij 9 x-^ mmjk ijk 

'*' g|^ ^mni jkl jkl ••• 

0 II _ II ^ 

^stmn “mnp ~ ^stmn "^mnp 


(57) 
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1 II 5^ I — I 

FT ^^stmn ^3cm®hn^fp3gq •^3x^3x~ 


+ _. 


9x^3Xg 


^chi jk' 


D.J.. .M.n [0]B.J 


+ [ J 


3Y^3x~ ^chijkl 

9^ n”.rmp.I.I. 9^ D„," jkl [01 ^ 


9^^ 0mnij[010ij + 

II 


ijkl 
II 


.] 


]} 


j. r 0 D — at F 

'"stmn ‘^mnpq ““^stmn ‘^mnpq 


(58) 


etc. 


in which the right hand side are determined by expanding the applied in a 
polynomial as in Eq. (48) and the equivalency equations used for any point 


P are 
I 


1 /:rP> 


II 


* I 


— P ^ 


^Oijkl ^^kl ) + ^kl (>" )> + e,^^ (r ) 

= ’ P 1'’ ^I 

II II _n - II _P 0 *II _n 

fzy\ . _ .r' - 


(59) 


^Oijkl (•” 0, , 


= -AC., 


II 


kl 


^-P 


'ijkl "-kl 


(r^) 


ijkl ^kl 


(r ) ; P in 


(60) 


The notation D-jjkl [01 > 0-j jkl [01 ••• '’'’^sn that the D's are evaluated 

at the point "0". The D's are defined in Eq. (39,40) ... . The super- 
scripts I and II are referred to the regions and Jlji occupied by the 
inhomogenei ti es . 
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INTERACTION ENERGY 

If an elastic body is subjected to surface tractions t^-(r,t) and 
body forces f^{7,t) , the induced elastic fields depend upon the elastic 
moduli of the body. Let (r,t) , e^-j (r,t) , u.j (r,t) and a^-j(r,t) , 

•(>">t) , u^- (r,t) be the stress, strain and displacement fields induced 

in _ f _ 

when the elastic moduli are C^-j|^-]{r) and C^-j|^-](r) , respectively. The 
initial state, denoted by superscript i£ , can be considered as the state 
where there are no inhomogeneities and the final state, denoted by f , 
can be considered as the state where there are inhomogeneities present. 

Using the notations given in [14], the difference in power input and 
the rate of change in kinetic energy plus potential energy is 


AE = - k""" + 


- J fi(ui^ - uj")dv - J t.j(u^-^ - u-")ds 
V s 


(61) 


where AE = interaction energy due to the presence of the inhomogeneities 

If * ■ 

K = J J p u,- u^- dv = kinetic energy 
^ v 

1 


U = 2 ! *^ij ^ij ~ strain energy 


Hence the interaction energy rate is 

•f *in f ,"f*f ..in»in, 

e' - E^" = / p(u| u.j^ - u. u. )dv 


1 ^ f 

V 


in • in. 


+ — f(a..u. . - a. . u. ,)dv 
2 V IJ ^ ,J TJ 




- t.(u/ - 


(62) 



21 


The equations of motion for the induced fields are. 

f .. f in " 

f in 

t. = 0 ..n. = 0.. n. 

1 U J J 

Using integration by parts and the equations of motion plus 
boundary conditions, Eqs. (63,64), the following identity can be 


derived: 


f'-f in'in.. 

■ // “/ - “ii" “,-'">"3 

- fff (a- ■ - a- U.''")dv 

'‘"iJ.J i iJ,3 1 

= ff t.(u.'*' - u-''")n ds 

s' -' ' J 

- m P(^r 
+ J/J 

V 

- |i K-r V">"o 

- W “° 13!3 " ■ ‘'’ 13.3 - 

- HI pC=," p," - "l’")'’' 

= HI (P,r “<(3 - “33- 

V 

+ HI - “V" 

- HI Cp u/ p/ - p-> p,’")-» 

V 

The substitution of Eq. (65) in Eq. (62) leads to 
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AE = - /J t.(u.^ - u.^")ds - 1 /// - u^^")dv 

S V 

1 1 iij pC«v 

■ - ill W 

- Iff Kf, - pii'/jp/ - (f, - pii,‘")Uj'"l dv 

- -jlll <“<]" “ns - “ilS)-*' 

1 rrr r f ' in in * f , . 

2 ij i,j IJ i.j 

The above equation may be expressed in terms of strains as follows: 

i i III (p,i p,; - p,i" p,5")d. 

■ -i III <“ij" si - si s-I")^» 


( 66 ) 


( 67 ) 


When inhomogeneities exist in the body, the following definitions are 
noted: 


• in *a in a _o *a 
'"ij " ^ij ’ °ij " '^ij “ ^ijkl ^kl ’ 


• f • a 

e . . = e . . + e . . 

iJ IJ ij 


(68) 


f 

7 . . 

TJ 


= 

^ijkl 

f a 
^^kl 


in matrix 


= 

1 jkl 

/ a 
^^kl 

^ "kl^ 

0 

" ^ijkl 

1 a 

'Si * Si 

* ^ 
- "kl ) 

= 

II 
^i jkl 

r ^ 

f^kl 


= c.°., 

ijkl 

'Si * Si 

- e* II 
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Hence the interaction energy due to the presence of inhomogeneities is 

AT r"E ri^ * 3. 

4E . E - E .--/// d» 

★ 

where n = + ... and C|^.| is the eigenstrain in each fl . 
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CONCLUDING REMARKS 

Quantitative determi nation of attenuation and velocity factors 
requires the solution of far fields and near fields for a material 
containing inhomogeneities under the excitation of incident power. The 
brief review of the existing scattering theory of a single flaw led to 
the realization of the importance of finding displacement and strain 
fields inside inhomogeneities. They are currently not available. Defini- 
tion and general equations for time-harmonic displacement and strain fields 
in a pair of interacting inhomogeneities are given. 

The interaction problems are presented via the dynamic eigenstrain 
concept. This approach converts the problem of dealing with inhomo- 
geneous boundary conditions to that of dealing with an nonhomogeneous 
differential equation. The nonhomogeneous term is directly related to 
the strains and displacements in the scatterer and can be obtained by 
the method of equivalent inclusions. 
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